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ABSTRACT 

Observations by the Cores to Disk Legacy Team with the Spitzer Space Telescope have identified a 
low luminosity, mid-infrared source within the dense core, Lynds 1014, which was previously thought 
to harbor no internal source. FoUowup near-infrared and submillimeter interferometric observations 
have confirmed the protostellar nature of this source by detecting scattered light from an outflow 
cavity and a weak molecular outflow. In this paper, we report the detection of cm continuum emission 
with the VLA. The emission is characterized by a quiescent, unresolved 90 /iJy 6 cm source within 
0'.'2 of the Spitzer source. The spectral index of the quiescent component is a = 0.37 ± 0.34 between 
6 cm and 3.6 cm. A factor of two increase in 6 cm emission was detected during one epoch and 
circular polarization was marginally detected at the 5a level with Stokes V/I = 48 ± 16% . We have 
searched for 22 GHz H2O maser emission toward L1014-IRS, but no masers were detected during 7 
epochs of observations between June 2004 and December 2006. L1014-IRS appears to be a low-mass, 
accreting protostar which exhibits cm emission from a thermal jet or a wind, with a variable non- 
thermal emission component. The quiescent cm radio emission is noticeably above the correlation 
of 3.6 cm and 6 cm luminosity versus bolometric luminosity, indicating more radio emission than 
expected. In this paper, we characterize the cm continuum emission in terms of observations of other 
low-mass protostars, including updated correlations of centimeter continuum emission with bolometric 
luminosity and outflow force, and discuss the implications of recent larger distance estimates on the 
physical attributes of the protostar and dense molecular core. 

Subject headings: radiation mechanisms:thermal, non-thermal — radio continuum:stars — stars: for- 
mation 
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1. INTRODUCTION 

It is extremely difficult to identify the incipient 
stages of low-mass (w 1 Mgun) star formation be- 
cause dense molecular cloud cores obscure nascent pro- 
tostars. Submillimeter dust continuum surveys (e.g., 
Ward-Thompson et al. 1994, Shirley et al. 2000, Visser 
et al. 2002, Kirk et al. 2005) have identified several 
dense cores with no apparent internal sources, based 
on the lack of an IRAS point source and the diffuse na- 
ture of submillimeter dust emission. Observations by the 
Cores to Disk Legacy Team (c2d) with the Spitzer Space 
Telescope have identified a few mid-infrared sources that 
are embedded near the submillimeter continuum peaks 
of previously classified starless cores (e.g. Young et al. 
2004, Bourkc et al. 2006). These new objects are of low- 
luminosity (Lint < 0.1 Lq) and presumably low-mass 
since they were not previously detected by IRAS. Some 
of these objects may be in the earliest stages of accre- 
tion. These newly identified low-mass, low-luminosity 
protostars warrant detailed follow-up studies to deter- 
mine their evolutionary status. 

The first newly identified object detected in the c2d 
survey, L1014-IRS (Young, et al. 2004), was mod- 
eled as a very low-luminosity (Lmt < O.ILq), low-mass 
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{M < O.IM©) object embedded within the Lynds 1014 
dark cloud (Lynds 1962) at a distance of approximately 
200 pc. This object has been classified as a VeLLO (Very 
Low-Luminosity Object) by the c2d team: an object with 
an internal protostellar luminosity < O.ILq that is di- 
rectly associated with a dense molecular core. The recent 
study of Morita et al. (2006) has suggested a revised dis- 
tance estimate of 400 to 900 pc based on the possible age 
ranges of nearby T-Tauri stars that are spatially within 
2° of the L1014 dense core. However, it is not clear that 
these T-Tauri stars are directly associated with L1014. 

Determining the evolutionary state of L1014-IRS has 
been the subject of several follow-up studies. The large 
scale molecular distribution in the dense core was deter- 
mined by the single-dish mapping survey of Crapsi et al. 
(2005). No evidence for a large scale CO outfiow was 
detected; however, a molecular outfiow was detected on 
smafl scales with the SMA (Bourke et al. 2005). Near- 
infrared observations detect scattered light, presumably 
from the outflow cone, at 1.6 /im and 2.2 /im (Huard et 
al. 2006). The SMA-detected CO outflow is ahgned with 
the direction of the near-infrared scattered light cavity. 
High resolution molecular observations with BIMA indi- 
cate that the protostar is not at the peak of the molec- 
ular and dust column density in the core, but offset by 
about 8" in the plane of the sky (Lai et al., in prepa- 
ration). This offset is also seen in (sub)millimeter con- 
tinuum maps (Young et al. 2004) and the near-infrared 
extinction map (Huard et al. 2006). 

Despite these significant observational efforts, a single, 
consistent picture of the evolutionary state of L1014-IRS 
has not emerged. In order to better characterize the 
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physical nature of L1014-IRS, we have conducted cen- 
timeter radio continuum observations using five array 
configurations of the Very Large Array ^ (VLA). Cen- 
timeter radio continuum emission is well correlated with 
the luminosity of protostellar sources (Anglada 1995) and 
is thought to arise from shock ionization from protostel- 
lar winds (Ghavamian & Hartigan 1998), from interac- 
tion of the protostellar jets with dense gas in the interface 
of the outflow cavity (Curiel et al. 1987, 1989, Shang et 
al. 2004), or from accretion shock-driven photoionization 
(Neufeld & Hollenbach 1996). 

In this paper we report the detection and characteriza- 
tion of the cm radio continuum emission toward L1G14- 
IRS (§3.1). We compare the detected cm emission with 
observations of other low-mass protostars (§4.1). We dis- 
cuss the non-detections of 22 GHz H2O masers and com- 
pare our upper limits to the recent maser monitoring 
surveys of low- mass protostars (§4.2). We also compare 
the derived source properties from the diverse studies of 
the protostar and dense core in terms of the range of 
distance estimates to L1014 (§4.3). 

2. VLA OBSERVATIONS 

L1014-IRS was observed during 10 epochs in 5 ar- 
ray configurations (D, A, BnA, B, and C) with the 
Very Large Array (Table 1). All observations were 
centered on the published Spitzer mid-infrared source 
{a = 2l''24™07^51, S = -F49°59'09'.'0, J2000.0). Con- 
tinuum observations were made at 3.6 cm, and 6.0 cm, 
with two polarization pairs at adjacent frequencies, pro- 
viding a total equivalent bandwidth of 172 MHz. We 
also attempted to detect H2O masers by observing the 
JxaKc = 61 6 52 3 transition at 22.23508 GHz with, 
typically, 24.4 kHz spectral resolution (0.3 km/s) span- 
ning ±20 km/s velocity coverage. 

For the 3.6 and 6.0 cm data, the data were reduced 
independently using the standard routines in AIPS+ + 
and AIPS. Complex gain calibration was performed by 
switching to the nearby quasar 2137+510, 2.4° from 
L1014-IRS, on time-scales of 15 to 30 minutes (Hamaker, 
Bregman, & Sault 1996a,b). The absolute fiux density 
and bandpass calibration were determined from observa- 
tions of the quasars 3C48 and 3C283. 

The Stokes / and V images were deconvolved using 
the Cotton-Schwab algorithm (e.g., Schwab 1984) and 
Clark-Hogbom algorithm (Hogbom 1974, Clark 1980) 
with a few thousand iterations and interactive CLEAN 
regions. Imaging the L1014-IRS field was difficult due to 
the presence of several bright sources within the VLA pri- 
mary beam (Figure la). Special care had to be taken in 
the CLEANing process (e.g., Cornwell, Braun, & Briggs 
1999) and multiple reductions with variations in the 
CLEAN parameters were performed. We have checked 
the consistency of our images by also reducing the data 
with the standard AIPS routines, and the fluxes agree 
within the statistical errorbars. Generally, the images are 
made with natural weighting of the visibilities (Briggs 
1995); however, uniform weighting was used to obtain 
better angular resolution for the full track (9 hour) ob- 
servations on the days of July 1, 2004 and August 21, 

^ The VLA is operated by NRAO. The National Radio Astron- 
omy Observatory is a facility of the National Science Foundation 
operated under a cooperative agreement by Associated Universi- 
ties, Inc. 



2004. 

3. RESULTS 

3.1. Radio Continuum Detections 

We detected cm radio continuum emission from a 
source within 0'.'2 of the Spitzer mid-IR source using the 
VLA at 3.6 cm and 6 cm (Figure 1). The initial detec- 
tions were made during 9 hour tracks at 3.6 cm and 6 
cm with the VLA in the D configuration. Subsequent 
observations detected the source at 6 cm in the other 
three VLA configurations, and again at 3.6 cm in D con- 
figuration. All of the 6 cm observations, except for the 
initial 6 cm detection on August 21, 2004, indicate a 
constant flux density source with an average 6 cm flux of 
88 ± 11 /iJy (80"). The source is unresolved in all array 
configurations. The two detections at 3.6 cm are also in 
agreement with each other despite being separated by 17 
months. The average 3.6 cm flux is 111 ± 8/uJy (14(j). 
The spectral index between two wavelengths (A2 > Ai) 
is defined as 



a ■ 



ln(A2/Ai) 



We find an average spectral index between 3.6 cm and 6 
cm of a = 0.37 ±0.34. 

We may estimate the likelihood that the radio contin- 
uum emission is directly related to L1014-IRS and not 
due to a background object by calculating the number 
of background extragalactic sources with sufficient flux 
that are expected within our synthesized beam. The cm 
continuum source is unresolved in the A configuration 
observations at 6 cm with a 0'.'46 x 0'.'44 synthesized 
beam (natural weighting). The position of the 6 cm 
source {a = 2l''24"07^53, S = +49°59'08'.'9, J2000.0) 
is within 0'.'2 of the Spitzer mid-infrared source. The 6 
cm source is the only source with a flux > 90 ^Jy within 
2' of the Spitzer source. The total number of background 
sources with a flux greater than 5* /iJy at 6 cm is given 



by 7V(> S) = 0.42(5'/30/iJy) 



-1.18 



per square arcmmute 



(Fomalont et al. 1991); therefore, the number of back- 
ground sources with a flux > 90 /zJy is 0.115 arcmin"^. 
The probability that a background source lies within our 
6 cm synthesized beam is 0.0006%. We believe the cen- 
timeter emission is associated with L1014-IRS and not 
associated with a background source due to the excep- 
tional positional coincidence and strength of the emis- 
sion. 

The cm continuum emission was not constant at all ob- 
served epochs. The 6 cm observations during the epoch 
of August 21, 2004 indicate a factor of 2 higher flux. 
The 6 cm detection is an llcr detection of 173 ± 16/uJy 
(Figure Ic). This indicates that the centimeter flux is 
variable and may be an indication of a flare (§4.1.2). We 
checked the consistency of the August 21 observations 
by comparing the 6 cm flux of other unresolved sources 
within the L1014 field during the other 6 cm epochs. No 
systematic increase was detected among the background 
unresolved sources on August 21, 2004. 

We checked for variability on short timescales by sep- 
arating and re-imaging the August 21, 2004 6 cm obser- 
vations into four 2 hour time blocks. The source was 
found to not vary, within the uncertainties, during each 
2 hour time block. This does not rule out a variation 
on shorter timescales; however, it is extremely difiicult 
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to detect that variation due to the higher noise level in 
progressively shorter time bloc;ks. 

A Stokes V source (circularly polarized) was 
marginally detected at the Scr level at 6 cm on August 
21, 2004 with a flux of 84 ± 17 /uJy at the position of the 
6 cm Stokes / source (Figure Id). The fraction of circular 
polarization, fc = Stokes V/I — 48 ± 16%, is quite high 
if the detection is significant. Since the Stokes V result 
was obtained during a full single track of VLA observa- 
tions, it will be difficult to confirm until L1014-IRS is 
observed with wider bandwidth (e.g., with the eVLA). If 
the circular polarization signal is real, then this observa- 
tion indicates that the radio emission observed on August 
21, 2004 must originate from a non-thermal mechanism 
(§4.1.2). 

3.2. Water Maser Search 

We searched for H2O masers during seven epochs span- 
ning 22 months by observing the Jk^k^ = Q16 523 
transition at 22.23508 GHz. No H2O masers were de- 
tected at any epoch. The combined Icr rms of the non- 
detection is 3.1 mJy/beam with a channel spacing of 24.4 
kHz and a total bandwidth of 3.125 MHz (~40 km/s). 
The individual observations are summarized in Table 1. 

4. DISCUSSION 
4.1. Centimeter Radio Continuum Emission 
4.1.1. Steady Component 

Centimeter continuum emission has been detected to- 
ward many but not all high-mass and low-mass proto- 
stars. The emission is most commonly thoiight to origi- 
nate from bremsstrahlung (free-free) emission from ion- 
ized gas, although some protostellar objects also display 
non-thermal emission. For high-mass protostars, the ion- 
ization mechanism is photoionization usually in the form 
of an embedded HII region (Churchwell 1990). For pro- 
tostars that are later in spectral type than B, the ioniz- 
ing radiation from the star is not enough to significantly 
photoionize the surrounding envelope and an alternative 
mechanism is needed to explain the observed emission 
(Rodriguez et al. 1989, Anglada 1995). Since nearly 
all low-mass, embedded protostars are known to have 
molecular outflows (e.g. Wu et al. 2004), the ioniza- 
tion is postulated to arise from shocks generated from 
a jet (e.g., Cohen, Bieging, & Schwartz 1982, Bieging 
& Cohen 1989, Curiel et al. 1987, 1989, Rodriguez & 
Reipurth 1996, Shang et al. 2004). Direct evidence for 
this hypothesis comes from observations using interfer- 
ometers at high angular resolution. Elongated centime- 
ter continuum emission regions are observed with the 
same orientation as the large-scale molecular outflow to- 
ward a few protostars (e.g., Anglada 1995; Bontemps, 
Ward-Thompson, & Andre 1996). In addition, the ob- 
served outflow force {Mq km/s/yr) theoretically provides 
enough energy in the shock to explain the observed cen- 
timeter fluxes toward most low-mass protostars (Cabrit 
& Bertout 1992, Skinner et al. 1993, Anglada 1995). 
Finally, the radio spectral index of many protostellar 
sources is consistent with optically thin (a ~ —0.1) to 
partially optically thick free-free emission between 3.6 
cm and 6.0 cm (2.0 > a > —0.1, e.g., Anglada et al. 
1998, Beltran et al. 2001). 

The centimeter continuum luminosity (e.g., Z/3.6 = 



S^.eD^ mjy kpc'^) of low-mass and intermediate-mass 
protostellar sources was first cataloged from the litera- 
ture in the review by Anglada (1995). Anglada plotted 
the 3.6 cm luminosity against the bolometric luminosity 
of protostars with Ltoi < 10"^ Lq and found a well corre- 
lated relationship (r = 0.79), L3.6 = 10-2-i(L6o(/lLo)°-'^ 
mJy kpc^. This relationship has formed the basis for 
predictions of the amount of centimeter emission that is 
expected in searches for new low-mass protostars (e.g., 
Harvey et al. 2002, Stamatellos et al. 2007). The 3.6 
cm luminosity correlation is directly related to the well 
established correlation of outflow force vs. bolometric lu- 
minosity (Bontemps et al. 1996, Wu et al. 2004); higher 
luminosity sources drive more powerful outflows that re- 
sult in a larger degree of shock ionization and therefore 
a larger 3.6 cm luminosity (Anglada 1995). 

Since 1995, many more centimeter observations have 
been made and the spectral coverage of the photometry 
of protostars has increased. We have used the detailed 
summary tables of Furuya et al. (2003; Table 4) and 
Anglada (1995) supplemented by the surveys of Eiroa 
et al. (2005) and Anglada et al. (1998) to catalog the 
3.6 cm, 6.0 cm, and bolometric luminosities of detected 
protostellar sources. We updated the bolometric lumi- 
nosity of sources observed in the submillimeter surveys 
of Shirley et al. (2000), Mueller et al. (2003), and Young 
et al (2003). The resulting sample of 58 sources at 3.6 
cm and 40 sources at 6.0 cm are plotted in Figure 2. We 
find updated correlations of 

logC-c-a.e/i "Jy ^pc^) = -(2.24 ± 0.03) + {0.71 ± 0.01) logct^^,/! Lq) (2) 

log(I,6.o/l mJy kpc^) = -(2.51 ± 0.03) + (0.87 ± 0.02) log(i:.j„,/l Lg) , (3) 

with correlation coefficients of r = 0.66 and r = 0.74 
respectively. This sample is not complete as there are 
many more protostellar centimeter detections for which 
no Lboi has been published. Nevertheless, we have up- 
dated the correlation of Anglada (1995) with twice as 
many points at 3.6 cm and plotted the correlation at 6.0 
cm for the first time. 

For comparison, wc have plotted the 3.6 cm and 6.0 cm 
luminosities of L1014-IRS in Figure 2 at the distances of 
200, 400, and 900 pc. L1014-IRS is above the correlation 
at all distances indicating that we have detected more 
centimeter continuum flux than expected. 

The spectral index of sources between 3.6 cm and 
6.0 cm is used to argue for the interpretation that the 
emission mechanism is consistent with partially optically 
thick free-free emission. Optically thin free-free emission 
is expected to have a spectral index a = —0.1 at cen- 
timeter wavelengths. In the optically thick limit, a ap- 
proaches 2.0 with intermediate values indicative of par- 
tially optically thick plasmas. We have also plotted the 
spectral index of sources from the literature that have 
been detected at both wavelengths and have a published 
Lboi determination in Figure 2. No correlation of a is 
observed with bolometric luminosity, probably indicat- 
ing that the optical depth associated with the jet's shock 
ionization is not dependent on the total protostellar lu- 
minosity or the strength of the molecular outflow. The 
median spectral index is a = 0.5, with most protostellar 
sources having flat or positive spectral indices. This me- 
dian value is close to the result expected for an ionized 
wind or jet with a 1/r^ densitv gradient (e.g., Panagia 
& Felh i975, Wright & Barlow 1975, Reynolds 1986). 
Unfortunately, most of the sources in Figure 2 were not 
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observed at both wavelengths on the same day and vari- 
abihty may result in significant scatter in the plot. 

The spectral index of the quiescent emission of L1014- 
IRS (0.37±0.34) is consistent with ionized free- free emis- 
sion with a density gradient. The spectral index agrees 
well with the median of the ensemble of protostellar 
sources measured. 

We have also updated the correlation between outflow 

force, Font (Mokm/s/yr), and centrimetric luminosity of 
Anglada (1995) using the molecular outflow compilations 
ofBontemptsetal. (1996), Furuya et al. (2003), and Wu 
ct al. (2004). The updated correlation of 44 sources is 
weak (r = 0.55), 

losC^out/l Mgkm/a/yr) = - (3.15±0.07) + (0.67±0.03) logCLg g/l mjy kpc^) . 

(4) 

This correlation is interpreted as evidence that jets from 

the molecular outflow provide enough shock ionization to 
account to the observed centimeter continuum emission 
(e.g., Anglada 1995). Assuming maximum ionization ef- 
ficiency, the minimum outflow force needed to create the 
observed level of 3.6 cm flux was estimated by Curiel et 
al. (1987, 1989) to be F^ut = 10~^-^L:i,e/l mJykpc^). 
This level is shown as a dashed line in Figure 2d. The 
observed outflow force toward L1014-IRS ranges from 
0.04 - 2.9 X 10-*^ Mq km/s/yr for distances of 200 to 
900 pc (see Bourke et al. 2005). At a distance of 200 pc, 
the upper limit on the outflow force is a factor of 2 lower 
than the Curiel theoretical minimum outflow force. The 
upper limit in the observed outflow force includes esti- 
mates of the missing flux due to interferometric spatial 
filtering as well as the average opacity corrections for 
low-mass protostellar outflows (see Bourke et al. 2005). 
Given the range of uncertainty in these estimates, the 
observed outflow force is below, but not necessarily in- 
consistent with the theoretical minimum outflow force 
needed to produce the observed 3.6 cm luminosity. The 
disagreement between the observed outflow force is more 
pronounced at larger distances, increasing to an order of 
magnitude below the theoretical curve for a distance of 
900 pc. 

The small observed molecular outflow force and the 
large observed centimeter continuum luminosities may 
indicate that another ionization mechanism is operating 
in L1014-IRS. There are a few other possibilities that 
have been discussed in the literature. We shall analyze 
the viability of two popular possibilities. 

Radiative transfer modeling of emission at 24 and 70 
fxm indicate a flux excess from a disk around L1014-IRS 
(Young et al. 2004). Neufeld & HoUenbach (1996) pos- 
tulated that the supersonic infall of material onto a pro- 
tostellar disk will create an accretion shock with enough 
ionization to generate ^ 1 mJy of continuum emission at 
centimeter wavelengths at distances of w 200 pc. How- 
ever, this mechanism does not appear to be able to pro- 
vide enough ionization to explain the emission observed 
toward L1014-IRS. In order to produce 90 /uJy emission 
at 6 cm at a distance of 200 pc, the protostellar mass has 
to be > 2 M0 and the accretion rate onto the disk must 
be > 10"'' Mo/yr (sec Figure 2 of Neufeld & HoUenbach 
1996). Estimates of the protostellar mass are very un- 
certain and highly distance dependent, however, 2 M0 is 
likely larger than the L1014-IRS protostar and disk mass, 
even for the far distance estimate of 900 pc (see Young 
et al. 2004). Furthermore, this accretion rate is an order 



of magnitude larger than the range of inferred accretion 
rates from the observed outflow momentum flux and the 
modeled internal luminosity of the source (< 3 x 10~^ 
M0/yr; Bourke et al. 2006, Young et al. 2004). Ion- 
ization from an accretion shock does not appear to be a 
likely explanation. 

A more plausible possibility is that there is a spheri- 
cal wind component. The expected continuum emission 
from self-shocked spherical winds have been modeled by 
numerous authors (e.g., Panagia & Felli 1975, Wright & 
Barlow 1975, Reynolds 1986, Gonzalez & Canto 2002). 
The recent study of Gonzalez & Canto model a time 
variable wind that generates internal shocks (Raga et al. 
1990) which produce ionization and centimeter contin- 
uum emission (see Ghavamian & Hartigan 1998). Their 
models produce centimeter continuum emission of sev- 
eral hundreds of /j,-Janskys and spectral indices that are 
positive for mass loss rates of 10~^ M0/yr at a distance 
of 150 pc. Accounting for the larger distance estimates 
of 200 to 900 pc and potentially lower mass loss rates for 
L1014-IRS, then this type of emission may still account 
for the flux observed at 3.6 cm and 6 cm (^ 100 /iJy). 
However, the spherical wind models is usually applied 
to evolved protostellar sources that are Class II (classi- 
cal T-Tauri stars) or later (e.g., Evans et al. 1987). If 
L1014-IRS is an older, more evolved protostar, then this 
may not be a problem (§4.3). 

4.1.2. Variable Component 

While the quiescent component of L1014-IRS is 
constant over 5 epochs with a positive spectral in- 
dex, the emission properties of L1014-IRS were signif- 
icantly different during the single epoch of August 21, 
2004. The 6 cm flux was larger by a factor of two, 
S'6.o(21AUG2004) = 173 ± 16 /xJy. Unfortunately, the 
spectral index of elevated emission was not determined 
since L1014-IRS was observed at a single wavelength. 
Circular polarization was detected at the 5(t level indi- 
cating non-thermal emission. Thus. L1014-IRS has vari- 
able centimeter emission, although the timescale of the 
variability is not constrained since it was seen to vary 
during only a single epoch. 

In general variability of centimeter continuum sources 
has not been properly addressed since observations of 
sources are limited to a few epochs. There has not been 
a systematic, monthly monitoring campaign of deeply 
embedded soiirces to characterize their centimeter vari- 
ability; however, there is observational evidence for vari- 
ability among deeply embedded protostellar sources. For 
instance, the Class source, B335, is known to vary be- 
tween an upper limit of < 80 /iJy (1994 December) and 
390 ^Jy (2001 January) at 3.6 cm (Avila et al. 2001, 
Reipurth et al. 2002). Another example is the variability 
and purported evidence for jet precession of the centime- 
ter continuum sources toward IRAS16293 (Chandler et 
al. 2005). Variability contributes scatter (de-correlation) 
of the luminosity correlations shown in Figure 2. In ad- 
ditional to the known variable thermal sources, several 
non-thermal protostellar sources are known to be highly 
variable (e.g. T Tauri stellar flares, see White 1996); 
but, most of those objects are more evolved than deeply 
embedded protostars. 

The observed increase in emission on August 21, 2004 
and 5cr Stokes V detection is indicative of variable non- 
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thermal emission toward L1014-IRS. While there have 
been a few high-mass protostars with observed nega- 
tive spectral indices (e.g., Reid et al. 1995; Garay et 
al. 1996), most embedded (< Class I) low-mass proto- 
stars with cm radio emission have positive spectral in- 
dices (see Figure 2). There are only a few embedded 
low-mass protostars toward which negative spectral in- 
dex, non-thermal emission is detected (e.g.. Shepherd 
& Kurtz 1999, Girart et al. 2002). One case, R CrA 
IRS 5, was detected with significant circular polarization 
(Feigelson et al. 1998). The authors postulate that the 
emission is due to gyrosynchrotron emission and may 
originate from magnetic reconnection events associated 
with flares (Feigelson et al. 1998). The physical mecha- 
nism for generating the radio flare is still not well under- 
stood (e.g., Basri 2004) and it is questionable whether 
it applic;s to the c;nibc;d(led phase of low-mass protostars 
(see below). A second case, IRAS 19243+2350, is a steep 
spectrum non-thermal source (a — —0.82 ± 0.04) that is 
elongated in the direction of its CO outflow (Girart et 
al. 2002). Girart et al. postulate that the emission may 
originate from a bi-conical synchrotron source tracing 
the protostellar jet, similar to observations of the high- 
mass som-ce W3(0H) (Reid et al. 1995, Wilner, Reid, & 
Menten 1999). Non-thermal emission may be present in 
low-mass protostellar jets, but the level of emission may 
be dominated by the thermal, shock-ionized component 
of the jet. In the case of L1014-IRS, this non-thermal 
jet component would have to be variable. We shall dis- 
cuss several possibilities for the origin of the observed 
non-thermal emission toward L1014-IRS. 

In order to better understand the non-thermal emis- 
sion mechanism, we estimate the brightness tempera- 
ture of the emission to be Tf, = S^^^cP /IkRl^^^^ = 
0.14 - 2.8 X 10^' K for distances of 200 to 900 pc and 
an emitting region that is i?cmit ~ 1 AU in size. The 
brightness temperature is very sensitive to the assumed 
size of the emitting region. A i?emit of 1 AU is appro- 
priate for a small flare; but, smaller when compared to 
the solar coronal emitting region which is typically less 
than K, 5 AU (e.g. Leto et al. 2000). Since the emission 
observed toward L1014-IRS was unresolved, even in the 
A-array configuration, then we can only limit the size of 
the emission region to < 90(-D/200pc) AU. For instance, 
if we assumed that the emitting region was 45 AU (half 
of our A-array resolution), then the brightness temper- 
ature drops to < 100 K. This is too low even for the 
steady, thermal free- free component (T ~ 10^ K), unless 
the emission was very optically thick. This cannot be 
the case since r >> 1 would imply a approaching 2.0 
which is not observed. While the size of the emitting re- 
gion is severely unconstrained, the detection of circular 
polarization indicates non-thermal emission probably on 
small (< few AU) size scales, most likely due to gyrosyn- 
chrotron emission (e.g., Ramaty 1969; Dulk & Marsh 
1982). 

Radio flares are observed toward very low-mass ob- 
jects including brown dwarfs (Berger et al. 2002, Os- 
ten et al. 2006), late M dwarfs (Berger 2006), and T- 
Tauri stars (Bieging & Cohen 1989; White, Pallavicini, 
& Kundu 1992). The typical level of quiescent emis- 
sion toward low-mass stars and brown dwarfs is 100 /zJy 
and the flares are 1 mJy (for nearby distances of « 30 



pc) with significant circular polarization (/c « 50% ) 
detected dining the fiaring events. The spectral lumi- 
nosity at 6 cm of L1014-IRS during the event is = 
AttD'^S^ = 4.7 X 10i5(£)/200pc)2 erg s'^ Hz"!. This 
is about 2400 times larger than the most luminous ob- 
served solar flares (Bastian 2004) and about 200 times 
larger than the flares observed by Berger (2006) toward 
late M dwarfs. If the radio emission is due to a flare, it 
must be a powerful flare since L1014-IRS is at least 20 
times farther away than the average distance of sources 
detected by Berger {{d) = 10.6 ± 5.1 pc). However, it is 
not larger than the typical non-thermal flaring emission 
observed toward yoimg T Tauri stars of L,y = 10^'"' — 10^^ 
erg s^i Hz^i (Giidcl 2002). The ratio of spectral flare 
luminosity to bolometric luminosity, L^/Li^oi = 4 x 10~^^ 
Hz~^, is also similar to the ratio observed toward classi- 
cal T-Tauri stars (see Giidel 2002). 

The timescale over which radio flares toward low-mass 
stars are observed tends to occur over minutes to hours. 
For instance, the low- mass brown dwarf, LP944-2, dis- 
covered by the 2001 NRAO summer students (Berger 
et al. 2002), displays flaring activity with an average 
timescale of 10 to 15 minutes. This is very different from 
the activity observed toward L1014-IRS on August 21, 
2004. We detected no evidence for short-term variabil- 
ity within our detected emission. The source appears to 
have a nearly constant flux that is twice as high as the 
steady component for at least an 8 hour period. This 
elevated emission then appears to be longer in duration 
than that observed during flaring events toward low-mass 
(proto)stars (Giidel 2002). 

An alternative possibility is that the elevated emission 
is not due to a flare, but due to rotational modulation of 
a non-thermal component associated with the magnetic 
connection between disk and accretion onto the star (i.e., 
Bieging & Cohen 1989). Such a mechanism has been 
postulated for the T-Tauri star, V410 Tauri, with a rota- 
tional period of 1.9 days. The observed emission toward 
V410 is 1 mJy with a negative spectral index. If the 
accretion spot is blocked from view for a fraction of the 
stellar rotation period, then it is possibly that we could 
have observed the source with the accretion spot is in 
view on August 21, 2004 and with the accretion spot 
blocked from view during the other epochs. The rota- 
tional period of L1014-IRS must be longer than 8 hours 
since elevated emission was observed during the entire 8 
hour track. A negative spectral index was observed to- 
ward V410 Tauri, while a negative spectral has not been 
observed toward L1014-IRS. This hypothesis is highly 
speculative and would require a regular monitoring cam- 
paign to test. 

Unfortunately, it is currently not possible to strongly 
constrain the origin of the non-thermal component to- 
ward L1014-IRS. The expanded bandwidth of the eVLA 
is needed to permit a systematic monitoring campaign 
with high enough signal-to-noise in only a few hour ob- 
servations to routinely check for a Stokes V detection 
and to determine an instantaneous spectral index. 

4.2. Water Maser Non-detections 

A compact, weak molecular outflow has been detected 
toward L1014-IRS (Bourke et al. 2005); therefore it may 
be possible to detect water masers if the jet is impinging 
on dense knots of material near the protostar. Previous 
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monitoring of water mascrs around low-luminosity pro- 
tostars (L < IOL0) indicate that the maser activity is 
highly variable and fugacious (e.g., Wilking et al 1994, 
Claussen et al. 1996, Furuya et al. 2003). There appears 
to be a threshold in luminosity {Liras > 25 Lq) above 
which water masers are always detected (Wilking et 
al. 1994). In general, low-luminosity protostars display 
less water mascr activity and have lower isotropic water 
maser luminosities than their high-luminosity counter- 
parts (e.g., Brand et al. 2004, Furuya et al. 2003). The 
lowest luminosity protostar with a water maser detection 
is the Class source, GF9-2, with Lboi = 0.3 L0 (Furuya 
et al. 2003). If the distance to L1014 is 200 pc, then 
protostellar luminosity of L1014-IRS is a factor of 3 less 
luminous than GF9-2. A search for water masers toward 
VeLLOs significantly expands the probed luminosity pa- 
rameter space. 

Water maser activity for low-mass protostars is quan- 
tified in terms of the isotropic maser luminosity, L^h°o — 
AttD"^ J S^dv (e.g. Wouterloot & Walmsley 1986). We 
can calculate the 3cr upper limit to the isotropic maser 
luminosity of L1014-IRS using the combined rms of 
the non-detections (3.1 mJy/beam in a 24.4 kHz chan- 
nel). The average linewidth of H2O masers sources with 
Lboi < 10 L0 is (Aw) = 1.2 ± 0.4 km/s (Furuya et al. 
2003), corresponding to 3.5 charmcls in the VLA spec- 
trum. Therefore, the 3ct upper limit to the maser in- 
tensity toward L1014-IRS, assuming the average maser 
linewidth, is 11.2 mJy/beam. 

Furuya et al. find a correlation between the average 
isotropic II2O mascr luminosity and Lhoi] however, there 
is a mistake in their conversion from integrated intensity 
(K km/s) to isotropic water maser luminosity (Lq) due 
to confusion of kpc with pc (Furuya 2007, private com- 
munication). The correct conversion equation is given 

by 



L5J°o = 2.5 X IO-^Lq 



(I 



Tldv 



V200pcy VlKkm/s 



(5) 



We have updated the Furuya et al. correlation by cor- 
recting the conversion factor in their Table 2 and also 
adding in points from the literature, namely from a sur- 
vey of high-mass star-forming regions and UCHII regions 
by Wouterloot & Walmsley (1986), Palla et al (1993), 
and a survey of Bok globules by Gomez et al. (2006). 
Unfortunately, we were not able to add points from some 
prominent water maser surveys such as Wilking et al. 
(1994) and Valderetto et al. (2002) due to incomplete re- 
porting or differences in the definition of the integrated 
flux. The updated correlation for 73 maser sources is 
shown in Figure 3. While at any given epoch, the wa- 
ter maser luminosity of a source is highly variable, the 
mean isotropic water maser luminosity is well-correlated 
with Lboi {r ~ 0.88) over 6 orders of magnitude in bolo- 
metric luminosity: LH2O = 3 x 1O~^L0 LjJ'^j'*. The re- 
lationship is less well correlated (r = 0.61), but docs 
not significantly change, if we restict the linear regres- 
sion to just low-luminosity sources with Lboi < 100 Lq, 
LH2O = 3 X 1O~^L0 L°oi^. The best fit slope is near the 
mean value from the many determinations made in the 
literature: e.g., Wouterloot & Walmsley (1986, L'jl.), 
FelU et al. (1992, L}P2), Palla et al. (1993, L?|^s), 



Brand et al. (2004, ). 



If we extrapolate this correlation to the internal lu- 
minosity of L1014-IRS at 200 pc (O.O9L0), then the ex- 
pected Lji^O is 3 X 10~^° L0. Therefore, our 3cr upper 
limits {Lh^o < 9.3 x 10-^'^{D/200pc)'^ Lg) are more 
than an order of magnitude below the predicted corre- 
lation and are significant compared to the strength of 
other detected low-luminosity protostars. The VLA is a 
very sensitive instrument for water mascr searches, and a 
systematic monitoring campaign toward newly identified 
VeLLOs with molecular outflows is needed to constrain 
the threshold for the excitation of water masers. 

4.3. The Distance and Evolutionary State of LIOI4-IRS 
4.3.1. Distance 

The evolutionary state of L1014-IRS depends on many 
observational factors and on the assumed distance to the 
protostar. Distances to dark clouds are very difficult to 
determine by standard star counting or extinction tech- 
niques (e.g., Schmidt 1975). For the previously published 
papers from the Cores to Disk Legacy group, a distance 
of 200 pc was assumed. Recent observation of T-Tauri 
stars, presumed to be associated with the L1014 dense 
core, indicate that the distance is greater than 400 pc 
and may be as large as 900 pc (Morita et al. 2006). The 
association of the T Tauri stars with the dense core is 
not clear as L1014 is located in a region with multiple 
dark clouds within a few degrees. Nevertheless, in light 
of these new results, we discuss the physical properties 
derived toward the L1014 core and L1014-IRS with dis- 
tances that range from 200 to 900 pc (see Table 2). 

One potential tool to constrain the distance to L1014 is 
the use of stellar color excess, E{B — V), of stars in the di- 
rection of L1014 (see Bourke et al. 1995). The MK Clas- 
sification Extension Catalogue (Morris-Kennedy 1983) is 
used to determine the spectral type and the intrinsic stel- 
lar color {B — V)q. If the dark clouds near L1014 are as- 
sociated with each other, then the extinction measured 
toward the MK stars may be used to look for discrete 
jumps in extinction with distance. This method is used 
to produce Figure 4 for stars within 10° of L1014. As 
evinced by E. E. Barnard's original photographs of this 
region (1927), there are many dark clouds within 10° of 
L1014. One of the clouds, B164 (3.8° from L1014), is es- 
timated to be 200 to 240 pc away based on its apparent 
association with the nearby star 80 Cygni (Pagani et al. 
1996). Another cloud, B361 (3.1° from L1014), was esti- 
mated to be 350 pc away (Schmidt 1975). Unfortunately, 
many of these clouds are probably not directly associated 
with L1014. As a result, there is no clear jump in extinc- 
tion with distance, but instead the jump in extinction 
occurs over a range of values from 200 pc to 500pc. 

Since many of the derived properties have different de- 
pendences on distance in Table 2, we can attempt to use 
them to limit the range of distances to the core. The 
mass of the dense core has been determined from dust 
continuum observations (Young et al. 2004, Md oc D^) 
and molecular line observations (Crapsi et al. 2005, 
Mmr OC D). The virial mass is equal to the dust mass 
(using the same radius to calculate both quantities) at 
a distance of 500 pc; however, if we allow for the typi- 
cal uncertainties in assumed dust opacity, then the dust 
mass is usually within a factor of 2 of the virial mass for 
star-forming cores (e.g., Shirley et al. 2002, 2003). Such 
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a range varies the distance from 250 pc to 1000 pc. 

A better constraint is found from the measured outer 
radius of the dense core. The radius was determined by 
near-infrared extinction observations of background stars 
(Huard et al. 2006). The size of low- mass star forming 
cores is hmited to typical outer radii of < 60, 000 AU 
(Shirley et al. 2002, Jorgensen et al. 2002, Kandori et 
al. 2005), or < 0.3 pc. The measured outer radius is less 
than 0.3 pc only for distances of 400 pc or less. Therefore, 
the size of the dense core argues for a closer value of the 
distance to L1014. 

The current set of observations are unable to place 
strong constraints on the distance to L1014, although 
there is evidence for a distance of < 500 pc. A detailed 
spectroscopic study of stars near the periphery of L1014 
is needed to better constrain the distance. 

4.3.2. Evolutionary State 

A fundamental question that remains to be answered 
is whether L1014-IRS is a newly formed protostar or a 
more evolved, accreting object. In order to address this 
question, we must synthesize the current observational 
information on L1014. The derived quantities are sum- 
marized in Table 2. 

First, there is a discrepancy between the modeled den- 
sity structure of the L1014 dense core from the near- 
infrared extinction map (Huard et al. 2006) and the 
(sub)millimetcr continuum maps (Young et al. 2004). 
The near-infrared data fit a Bonnor-Ebert sphere (Ebert 
1955, Bonnor 1956), a hydrostatic pressure-bounded con- 
figuration, with a central density that is a factor of 8 
higher than the (sub)-millimcter models. Both methods 
have biases that compromise the estimates of the cen- 
tral density. The (sub)-millimeter observations detect 
the core at less than lOcr, resulting in a relatively low 
signal-to-noise radial profile. The ability to successfully 
model the size of the "plateau" region of a Bonnor-Ebert 
profile is marginalized by the poor signal-to-noise and 
telescope resolution (Shirley & Jorgensen 2007, in prepa- 
ration). Therefore, the (sub)-millimeter models are a 
lower limit to the central density of the core. The near- 
infrared extinction maps arc limited by the number of 
background stars that can be traced at high extinction 
{Ay > 30 mag). In the case of L1014, there are only two 
stars above 35 mag that heavily bias the smoothed den- 
sity structure (T. Huard 2006, private communication). 
Therefore, near-infrared extinction mapping is unable to 
reliably determine the size of the "plateau" region for 
densities above a few 10^ cm~'^. The true central den- 
sity of L1014 probably lies between 10^ cm~'^ and 10^ 
cm~^. These value indicate a core that is moderately to 
strongly centrally condensed, but there are dense star- 
less cores that have higher central densities (e.g., L1544 
Evans et al. 2001; L183, Pagani et al. 2003) 

The evolutionary state of L1014-IRS may not be di- 
rectly tied to the density structure of the core since 
L1014 may have fragmented in the past. There is evi- 
dence for this hypothesis since L1014-IRS is not located 
at the peak of the dust column density or molecular col- 
umn density distributions (Young et al, 2004, Crapsi et 
al. 2005, Huard et al. 2006, Lai et al. 2007 in prepa- 
ration), unlike the case for other Class protostars (see 
Shirley et al. 2000, Jorgensen et al. 2002). The proto- 
star is located 8" to the north of the core peak or 1600 



to 7200 AU in separation in the plane of the sky for dis- 
tances of 200 to 900 pc. While this distance is within 
the core's FWHM contour (see Jorgensen et al. 2006), it 
may be evidence of fragmentation. Fragmentation of the 
core is a more likely possibility due to the higher than 
average (for low-mass, dense cores) level of turbulence 
observed toward L1014 (MacLow & Klesscn 2004); it is 
higher than observed in most starless cores with a CS 
hnewidth of 0.68 km/s (Crapsi et al. 2005; Lee, Myers, 
& Plume 2004). The protostar may have also moved from 
its birth-site near the peak of the core. This can occur if 
non-spherical accretion is considered (Stamatellos et al. 
2005). For a reasonable protostellar velocity of 0.1 km/s 
(Walsh et al. 2007), it would take L1014-IRS 75,000 to 
350, 000 years to move 8" in projection for distances of 
200 to 900 pc. These numbers are within the estimated 
lifetimes of the prestellar and Class protostellar phases 
(Kirk et al. 2005, Ward-Thompson et al. 2007). 

L1014-IRS has a disk (Young et al. 2004) and is 
currently accreting material at a relatively low rate 
based on its weak molecular outflow and low luminosity. 
The oiitflow has been active for at least several thou- 
sand years since an observable outflow cavity has been 
cleared over a region greater than 10" (corresponding to 
> 2000(£>/200pc) AU) and has a wide opening angle 
{6 > 100°, Huard et al. 2006). L1014-IRS is still classi- 
fied as a Class source; however, if the distance is greater 
than 200 pc, L1014-IRS would not be strictly classified 
as a VeLLO since its internal luminosity is most likely 
greater than 0.1 Lq. 

The radio continuum observations indicate that L1014- 
IRS is one of the lowest luminosity protostars toward 
which centimeter continuum emission is detected. The 
3.6 cm and 6 cm emission is higher than expected for 
its protostellar luminosity and the outflow force is lower 
than theoretically expected for its centimeter continuum 
luminosity. Shock ionization from the protostellar out- 
flow is marginally a plausible explanation for nearby dis- 
tances (D < 500 pc). Another possible explanation may 
be a self-shocked spherical wind, although this type of 
emission is usually associated with more evolved proto- 
stellar objects. Variable emission is observed and the 
properties are consistent with the flaring properties ob- 
served toward T-Tauri stars except for the duration of 
the elevated emission. The true nature of the variable 
non-thermal component is still a mystery. The charac- 
teristics of both the thermal steady component and the 
variable non-thermal component are consistent with cen- 
timeter continuum emission from the later stages of pro- 
tostellar evolution. The bulk evidence seems to indicate 
that L1014-IRS is not an extremely young protostar (e.g. 
a flrst or second hydrostatic core), but instead seems to 
have been accreting material for at least several thousand 
years to tens of thousands of years and is currently ob- 
served in a low accretion state because it has decoupled 
from the peak of the L1014 dense core. 

Future observations are needed to better character- 
ize the centimeter continuum emission variability and to 
better constrain the mass loss properties of L1014-IRS. 
Observations with the new eVLA correlator will permit 
"instantaneous" spectral index studies due to the greatly 
enhanced bandwidth. Deep monitoring (rms < 20 /uJy) 
of the centimeter emission of low-mass protostars will 
be possible with integration times less than 1 hour. If 
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L1014-IRS is observed during a "flaring" event with the 
eVLA, it will be much easier to measure the degree of cir- 
cular polarization and to test for short term variability. 
Near-infrared observations of the protostar with large 
aperture telescopes (e.g., Gemini 8-m) of accretion di- 
agnostics, such as Br7, are needed to better constrain 
the accretion properties of L1014-IRS. Finally, high an- 
gular resolution observations of the dense core with an 
interferometer with a wide range of m, v coverage and 
brightness sensitivity (e.g., eVLA, CARMA, ALMA) will 
permit detailed analysis of the density structure of the 
inner core and the properties of the core in the offset re- 
gions between the protostar and the core column density 
peak. 

5. CONCLUSIONS 

We have detected cm radio continuum emission from 
L1014-IRS at 3.6 cm and 6.0 cm. The emission is charac- 
terized by a quiescent, nearly constant unresolved com- 
ponent of 90 fiJy at 6 cm with a spectral index of 
a = 0.37 ± 0.34 between 3.6 cm and 6 cm. We have up- 
dated, using recently published observations of low-mass 
protostars, the correlations of 3.6 cm and 6.0 cm con- 
tinuum luminosity versus bolometric luminosity and the 
correlation of molecular outflow force with 3.6 cm lumi- 
nosity. The quiescent emission component is consistent 
with partially ionized free- free emission, but is above the 
linear correlation for low-mass protostars between 3.6 cm 
luminosity and bolometric luminosity and below the cor- 
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relation of outflow force and 3.6 cm luminosity. The qui- 
escent continuum emission may be explained by shock- 
ionization from the protostellar jet or a time-variable 
wind. A non-thermal brightening by a factor of 2 of the 
6 cm continuum was detected on August 21, 2004 with 
48 ± 16% circular polarization at the 5a level. The true 
nature of the variable non-thermal component remains 
a mystery although the ratio of the flare luminosity to 
bolometric luminosity is consistent with those observed 
toward T-Tauri flares. We have also updated the corre- 
lation of isotropic water maser luminosity with bolomet- 
ric luminosity. We do not detect water masers toward 
L1014-IRS, consistent with the highly variable nature of 
water masers around low-luminosity protostars. Analysis 
of the derived properties of the L1014 dense core and the 
L1014-IRS protostar indicate that it is probably not an 
extremely young protostar but is a low-luminosity source 
that appears to have been accreting for at least several 
thousand years and is currently in a low accretion state. 
Properties of the protostar and dense core indicate it is 
probably at a distance of < 500 pc. 
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18 


/iJy 


11 Jan 2005 


BnA 


22.23508 




0'.'12 X O'.'IO 


24.4 kHz 




9.5 


mJy 


07 Mar 2005 


B 


22.23508 




0'.'35 X 0'.'30 


97.6 kHz 




6 


mJy 


09 Apr 2005 


B 


4.86 




1'.'7 X 1'.'4 


172 MHz 


87 


15 


/iJy 




B 


22.23508 




0'.'37 X 0'.'30 


24.4 kHz 




6 


mJy 


19 Jul 2005 


C 


4.86 




5'.'0 X 4"4 


172 MHz 


89 


21 


MJy 


05 Dec 2005 


D 


8.46 




9('5 X 5" 9 


172 MHz 


119 


20 


MJy 


29 Nov 2006 


C 


22.23508 




1'.'15 X 0^97 


48.8 kHz 




4.9 


mJy 


28 Dec 2006 


C 


22.23508 




1'.'64 X 0'.'83 


48.8 kHz 




9.8 


mJy 


Combined 


ABC 


4.86 




I'.'e X 1'.'5 


172 MHz 


84 


10 


MJy 



All epochs are imaged with natural weighting except for 01 Jul, 2004 and 21 Aug, 2004 where uniform weighting was used. 
^ The total bandwidth for continuum observations and the channel spacing for H2O maser observations. 



Kandori, R., Nakajima, Y., Tamura, M., ct al. 2005, AJ, 130, 2166 
Kirk, J. M., Ward-Thompson, D., & Andre, P. 2005, MNRAS, 360, 
1506 

Leto, G., Pagano, I., Linsky, J. L., Rodono, & Umana, G. 2000, 

A&A, 359, 1035 
Lee, C. W., Myers, P. C, & Plume, R. 2004, ApJS, 153, 523 
Lynds, B. T. 1962, ApJS, 7, 1 

Mac Low, M.-M. & Klessen, R. S. 2004, Rev. of Mod. Phys., 76, 
125 

Morris-Kennedy, P. 1983, MK Classification Extension. Mt. 

Stromlo Obs., Canberra (also in The Astronomical Data Centre 

CD-ROM Selected Astronomical Catalogs, Vol. 1, 1989) 
Morita, A., Watanabe, M., Sugitani, K., Itoli. Y.. l'(-liara, M., 

Nagashima, C, Ebizuka, N., Hasegawa, T., Kiuugasa, K., & 

Tamura, M. 2006, PASJ, 58, 41 
Neufeld, D. A., & HoUenbach, D. J. 1996, ApJ, 471, L45 
Osten, R. A. & Jayawardhana, R. 2006, ApJ, 644, 67 
Pagani, L. & Breart de Boisanger, C. 1996, A&A, 312, 989 
Pagani, L., Lagache, G., Bacmami, A. 2003, A&A, 406, 59 
Palla, F., Cesaroni, R., Brand, J., Caselli, P., Comoretto, G., & 

FeUi, M. 1993, A&A, 280, 599 
Panagia, N. & FeUi, M. 1975, A&A, 39, 1 

Raga, A. C, Canto, J., Binette, L. & Calvert, N. 1990, ApJ, 364, 
601 

Ramaty, R. 1969, ApJ, 158, 753 

Reid, M. J., Argon, A. L., Masson, C. R., Menten, K. M., & Moran, 

J. M. 1995, ApJ, 443, 238 
Reipurth, B., Rodriguez, L. F., Anglada, G., & Bally, J. 2002, AJ, 

124, 1045 
Reynolds, S. P. 1986, ApJ, 304, 713 

Rodriguez, L. F., Myers, P. C, Cruz-Gonzalez, I.. & Terebey, S. 

1989, ApJ, 347, 461 
Rodriguez, L. F., & Reipurth B. 1996, RMxA&A, 32, 27 
Schmidt, E. G. 1975, MNRAS, 172, 401 
Schwab, F. R. 1984, AJ, 89, 1076 

Shang, H., Lizano, S., Glassgold, A., & Shu F. 2004, ApJ, L69 
Shepherd, D. S. & Kurtz, S. E. 1999, ApJ, 523, 690 
Shirley, Y. L., Evans, N. J., II, & Rawlings, J. M. C. 2000, ApJS, 
131, 249 



Shirley, Y. L., Evans, N. J., II, & Rawlings, J. M. C. 2002, ApJ, 

575, 337 

Shirley, Y. L., Evans, N. J., II, Young, K. E., Knez, C. & Jaffe , D. 

T. 2003, ApJS, 149, 375 
Skinner, S. L., Brown, A., & Stewart, R. T. 1993, ApJS, 87, 217 
Stamatellos, D., Whitworth, A. P., Boyd, D. F. A., & Goodwin, S. 

P. 2005, A&A, 439, 159 
Stamatellos, D., Ward-Thompson, D., Whitworth, A. P., & 

Bontemps, S. 2007, A&A, 462, 677 
Valdettaro, R., Palla, F., Brand, J., Cesaroni, R., Comoretto, G., 

FeUi, M., & Palagi, F. 2002. A&A, 383, 244 
Villuendas, E., Anglada, G., Estalella, R,., Rodriguez, L. F., Curiel, 

S., & Torrelles, J. M. 1996, in Radio Emission from the Stars 

and Sun, ASP Conf. Ser. 93, eds. A. R. Taylor & J. M. Paredes, 

56 

Visser, A. E., Richer, J. S., & Chandler C. J. 2002, AJ, 124, 2756 
Walsh, A. J., Myers, P. C, DiFrancesco, J., Mohanty, S., Bourke, 

T. L., Gutermuth, R., & Wilner, D. 2007, ApJ, 655, 958 
Ward-Thompson, D., Scott, P. F., Hills, R. E., & Andre, P. 1994, 

MNRAS, 268, 276 
Ward-Thompson, D., Andre, P., Crutcher, R., Johnstone, D., 

Onishi, T., & Wilson C. 2007, in Protostars and Planets V, 

eds. B. Reipurth, D. Jewitt, & K. Keil, University of Arizona 

Press, Tucson, 33 
White, S. M., Pallavicini, & Kundu, M. R. 1992, A&A, 259, 149 
White, S. M. 1996, in Cool Stars, Stellar Systems, and the Sun, 9th 

Cambridge Workshop, APS Conf. Ser., 109, eds., R. Pallavicini 

& A. K. Duprec, 21 
Wilking, B. A., Claussen, M. J., Benson, P. J., Myers, P. C, 

Terebey, S., & Wootten, H. A. 1994, ApJ, 431, L119 
Wilner, D. J., Reid, M. J., & Menten, K. M. 1999, ApJ, 513, 775 
Wouterloot, J. G. A. & Walmsley, C. M. 1986, A&A, 168, 237 
Wright, A. E., & Barlow, M. J. 1975, MNRAS, 170, 41 
Wu, Y., Wei, Y, Zhao, M. Shi, Y., Yu, W., Qin, S., & Huang, M. 

2004, A&A, 426 503 
Young, C. H., Shirley, Y. L., Evans, N. J., II, Rawlings, J. M. C. 

2003, ApJS, 145, 111 
Young, C. H., Jorgensen, J. K., Shirley, Y. L., et al. 2004, ApJS, 

154, 396 



10 



TABLE 2 

L1014 Derived Properties with a Distance Dependence 



Property 


Dependence 


D=200pc 


D=400pc 


D=900pc 


units 


Ref. 




D 


> 1600 


> 3200 


> 7200 


AU 


1,3,5 


Lint 




0.09 


0.36 


1.8 


L© 


1 




£|2 


1.7 


6.8 


34.4 


M© 


1 


c 


£1-1/2 


1.5 xlO^ 


1.1 xlO^ 


7.1 xlO^ 




1 




D 


4.2 


8.4 


19.0 


M© 


2 




D 


0.13 - 0.16 


0.26 - 0.32 


0.59 - 0.72 


pc 


3 


^out 




> 60° 


> 60° 


> 60° 




3 


^out 


D 


540 


1080 


2430 


AU 


4 






0.14 - 1.4 xlO-'' 


0.56 - 5.6 xlO-* 


0.28 - 2.8 xlO" 


-3 M© 


4 


.dyn 
out 


D 


700 


1400 


3150 


yr 


4 


Mout 


D 


< 2 X lO"''' 


< 4 x 10-'^ 


< 9 X 10-'^ 


Me/yr 


4 


P out 




0.24 - 4.5 xlO"* 


0.1 - 1.8 xlO"-'* 


0.5 - 9.1 xlO- 


3 M© km/s 


4 




132 


0.02 - 1.4 xlO-3 


0.1 - 5.8 xlQ--^ 


0.05 - 2.9 xlO- 


-2 M© km2 / s2 


4 


T mech 
^out 


D 


0.03 - 2.5 xlQ-* 


0.7 - 5.0 xlQ-" 


0.15 - 1.1 xlO" 


L© 


4 


-pobs 
out 


D 


0.34 - 6.4 xlQ-^ 


0.07 - 1.3 xlO"'' 


0.15 - 2.9 xlO" 


M© km/s/yr 


4 


L3.6 


D'' 


4.4 xlO"^ 


1.8 X 10"^ 


9.0 X 10-^ 


mjy kpc'^ 


6 


Ls.o 


132 


3.5 xlO-3 


1.4 X 10-2 


7.1 X 10-2 


mjy kpc2 


6 


J flare. 
'-•1/ 


i>2 


3.9 xlO^^ 


1.6 xlO^^ 


7.9 xlO^^ 


erg s-i Hz-i 


6 


T ISO 




< 9 X 10-12 


< 4 X 10-" 


< 2 X 10-1° 


L© 


6 


References. 


— 1. Young et al. 2004; 2. Crapsi et al. 2005; 3. Huard et al. 


2006; 4. Bourke et al. 


2006; 5. Lai et al., in 


preparation. 



6. This paper. 



Note. — Source properties: A(?coVe = planc-of-sky distance between L1014-IRS and dense core peak, l^int = L1014-IRS internal 
luminosity (47r_D2 J S*dv), Men,, = mass of the envelope, nJT^™ = Bonnor-Ebert central density from mm continuum, M„ir = envelope 

virial mass, R""" = outer radius. Rout = outflow radius, iout = outflow inclination, Mout = outflow mass, t^^" = outflow dynamical time, 
Mout = outflow mass-loss rate. Pout = outflow momentum, Eout = outflow energy, LJJIfj'^'' = outflow mechanical luminosity, F'^^ = outflow 
force. La g = 3.6 cm quiescent luminosity (^.g-D^), Lg = 6 cm quiescent luminosity {Sq.oD'^\ hl^"^" = 6 cm flare spectral luminosity 
(i-irD^Se.o), Lif°o ~ isotropic H2O maser luminosity [iirD'^ J S^dy). 
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Fig. 1. — The initial detection of L1014-IRS at 3.6 cm (upper left), ABC-array combined 6 cm image (upper right), 6 cm "flare" on 21 
August, 2004 (lower left), and Stokes V6 cm detection for 21 August, 2004 (lower right). The synthesized beam is shown in the lower left 
of panels (b) and (c). The 2a contours in each panel correspond to; (panel b) ±20 /iJy/beam; (panel c) ±33 /iJy/beam; (panel d) ±33 
/xJy/beam. Negative contours are shown as light, dash-lines. 




Fig. 2. — Top: Updated correlation of 3.6 cm and 6.0 cm luminosity vs. bolomctric luminosity for protostars. The solid line is a linear 
regression for all points (excluding L1014). The L1014-IRS 3.6 cm and 6.0 cm quiescent luminosities are shown as trangles for distances 
of 200, 400, and 900 pc (always increasing left to right in all panels). Bottom Left: The spectral index of sources from the literature with 
L1014^IRS plotted as trangles at different distances. The dashed red line is the <x expected for optically thin free-free emission. Bottom 
Right: The outflow force (Font = Pout/tdyn) plotted vs. centimeter continuum luminosity. The solid line is a linear regression to the points 
(L1014-IRS excluded) while the dotted line is the Anglada correlation from 1995. The dashed line is the theoretical minimum relationship 
from Curiel et al. (1989). L1014-IRS is plotted as triangles for the three standard distances. The errorbars in Font represent the limits 
determined by Bourke et al. (2005) and are not statistical errorbars. 
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w A palla et al. 1993 



□ Wouterloot et al. 1986 

y 

2 I I I I I I I I I I I I I I I I 

-2 2 4 6 

log kol (Lq) 

Fig. 3. — Updated correlation of isotropic water maser luminosity vs. the bolometric luminosity. The filled points are from Furuya et 
al. (2003), the open circle points arc from Gomez et al. (2006), the open triangle points arc from Palla et al. (1993), and the open square 
points are from Wouterloot et al. (1986). The circular filled points are sources that were detected at all epochs in the Furuya et al. survey. 
The square filled points had at least one non-detection while the triangle filled points only had one detection. The errorbars on the Puruya 
et al. points correspond to the maximal and minimal isotropic water maser luminosity observed. The solid line is the linear regression for 
all the points in the figure (excluding L1014-IRS). The upper limits in the lower left axe for VLA observations of L1014 at distances of 200, 
400, and 900 pc. 
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Colour Excess vs Distance 
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Fig. 4. — The color excess vs. distance of MK Extension Catalog stars within 10° of L1014. The distance modulus is derived assuming 
a ratio of selective to total extinction of Ry = 3.1 = Ay/E{B — V) (Weingartner & Draine 2001). 



